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ABSTRACT 

We present the rest-frame Js-band and Ks-band luminosity function (LF) of a sample 
of about 300 galaxies selected in the HDF-S at Ks^ 23 (Vega). We use calibrated 
photometric redshift together with spectroscopic redshift for 25% of the sample. The 
accuracy reached in the photometric redshift estimate is 0.06 (rms) and the fraction 
of outliers is 1%. We find that the rest-frame Js-band luminosities obtained by ex¬ 
trapolating the observed Js-band photometry are consistent with those obtained by 
extrapolating the photometry in the redder H and Ks bands closer to the rest-frame Js, 
at least up to z ~ 2. Moreover, we find no significant differences among the luminosi¬ 
ties obtained with different spectral libraries. Thus, our LF estimate is not dependent 
either on the extrapolation made on the best-fitting template or on the library of mod¬ 
els used to fit the photometry. The selected sample has allowed to probe the evolution 
of the LF in the three redshift bins [0;0.8), [0.8;1.9) and [1.9;4) centered at the median 
redshift Zm — [0.6,1.2, 3] and to probe the LF at Zm — 0.6 down to the unprecedented 
faint luminosities Mjg ~ —13 and ~ —14. We find hints of a raise of the faint 
end (Mjs > —17 and M/yg > —18) near-IR LF at Zm ~ 0.6, raise which cannot be 
probed at higher redshift with our sample. The values of a we estimate are consistent 
with the local value and do not show any trend with redshift. We do not see evidence 
of evolution from z = 0 to Zm ~ 0.6 suggesting that the population of local bright 
galaxies was already formed at z < 0.8. On the contrary, we clearly detect an evolution 
of the LF to Zm ~ 1.2 characterized by a brightening of M* and by a decline of (p*. To 
Zm ~ 1.2 M* brightens by about 0.4-0.6 mag and p* decreases by a factor 2-3. This 
trend persists, even if at a less extent, down to Zm ~ 3 both in the Js-band and in 
the Ks-band LF. The decline of the number density of bright galaxies seen at z > 0.8 
suggests that a significant fraction of them increases their stellar mass atl<z<2 — 3 
and that they underwent a strong evolution in this redshift range. On the other hand, 
this implies also that a significant fraction of local bright/massive galaxies was already 
in place at z > 3. Thus, our results suggest that the assembly of massive galaxies is 
spread over a large redshift range and that the increase of their stellar mass has been 
very efficient also at very high redshift at least for a fraction of them. 

Key words: Galaxies: evolution - Galaxies: luminosity function - Galaxies: forma¬ 
tion. 


1 INTRODUCTION 

The luminosity function (LF) of galaxies is a fundamen¬ 
tal statistical tool to study the populations of galaxies. Its 
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dependence on morphological type, on wavelength and on 
look-back time provides constraints on the evolution of the 
properties of the whole population of galaxies, of the popu¬ 
lations of various morphological types and on their contribu¬ 
tion to the luminosity density at different wavelength. The 
parameters derived by the best-fitting of the LF, the char- 
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acteristic luminosity M*, the slope a and the normalization 
(p, provide strong constraints to the models of galaxy for¬ 
mation. 

The recent estimates of the LF based on local wide sur¬ 
veys in the optical, such as the Two Degree Field (2dF; 
Folkes et al. 1999; Madgwick et al. 2002; Norberg et al. 
2002) and the Sloan Digital Sky Survey (SDSS; Blanton et 
al. 2003) and in the near-IR, such as the Two Micron All 
Sky Survey (2MASS; Kochanek et al. 2001), have provided a 
comprehensive view of the local LF for different morpholog¬ 
ical types and wavelengths. It is now well established that 
the LF depends on morphological type, that the faint-end is 
increasingly dominated by galaxies with late-type morphol¬ 
ogy and spectra and that their number density increases 
going to lower luminosities. On the contrary, the bright-end 
is dominated by early-type galaxies whose fraction increases 
with luminosities (e.g. Marzke et al. 1994, 1998; Folkes et 
al. 1999; Kochanek et al. 2001; Bell et al. 2003). 

The studies of the LF based on optically selected red- 
shift surveys at a < 1 have provided the hrst evidences of 
a differential evolution of galaxies. Lilly et al. (1995a), us¬ 
ing the Canada France Redshift Survey (CFRS; Lilly et al. 
1995b) show that the rest-frame B-band LF of the blue pop¬ 
ulation of galaxies brightens by about 1 mag to 2 ~ 1 con¬ 
trary to the red population which shows a little evolution 
over the redshift range probed. They used the observed I- 
band photometry to derive the rest-frame B-band luminosi¬ 
ties. Subsequent studies conhrmed the different behaviors 
followed by the LF of the different populations of galaxies 
at various redshifts (e.g. Lin et al. 1997; Liu et al. 1998) 
and the differential evolution they underwent (e.g. Bell et 
al. 2004; Wolf et al. 2003). At higher redshift {z > 1) the 
studies of the LF in the optical rest-frame confirm the pres¬ 
ence of the bi-modality and provide evidences of luminosity 
and of density evolution down to 2 ~ 2 — 3 (e.g. Poli et al. 
2003; Gabasch et al.2004; Giallongo et al. 2005). 

Contrary to the light at UV and at optical wavelengths, 
the near-IR light is less affected by dust extinction and by 
ongoing star formation e.g. Rix & Riecke 1993; Kauffmann & 
Chariot 1998). Moreover, since the near-IR light of a galaxy 
is dominated by the evolved population of stars, the near- 
IR light is weakly dependent on galaxy type and it is more 
related to the stellar mass of the galaxy. Therefore, the evo¬ 
lution of the LF of galaxies in the near-IR rest-frame can 
provide important clues on the history of the stellar mass 
assembly in the Universe rather than on the evolution of the 
star formation. 

The studies of the K-band luminosity function con¬ 
ducted so far have shown no or little evolution of the pop¬ 
ulation of galaxies at 2 < 0.4 — 0.5 (e.g. Glazebrook et al. 
1995; Feulner et al. 2003; Pozzetti et al. 2003) with respect 
to the local population (Kochaneck et al. 2001; Cole et al. 
2001; Glazebrook et al. 1995; Gardner et al. 1997). On the 
contrary, evidences of evolution emerge at 2 > 0.5, even if 
some discrepancies are present among the results obtained 
by the various authors (see e.g. Cowie et al. 1996; Pozzetti et 
al. 2003; Drory et al. 2003; Caputi et al. 2004, 2005; Dahlen 
et al. 2005). One of the possible reasons of these discrepan¬ 
cies is related to the rest-frame K-band luminosities, usually 
extrapolated from the observed K-band photometry making 
use of the best-fitting template which should reproduce the 
unknown SED of the galaxy. To minimize the uncertain¬ 


ties, photometry at wavelengths long-wards the rest-frame 
K-band would be required. Given the obvious difficulties in 
getting deep mid-IR observations, some authors have used 
the observed K-band magnitudes to construct the rest-frame 
J-band LF (e.g. Pozzetti et al. 2003; Bolzonella et al. 2002; 
Dahlen et al. 2005) and to constrain the evolution of the LF 
in the near-IR. Another reason could be that deep K-band 
observations allowing a good sampling of the LF at mag¬ 
nitudes fainter than M* over a wide redshift range are not 
so common. Consequently the LF faint-wards of the knee 
is progressively less constrained with increasing redshift af¬ 
fecting the LF estimate. 

In this work, we aim at measuring the rest-frame Js- 
band and Ks-band LF of galaxies and its evolution to 2 ~ 3 
through a sample of about 300 galaxies selected at Ks^ 23 
on the Hubble Deep Field-South (HDF-S). The extremely 
deep near-IR observations collected on this field coupled 
with the optical HST data have allowed us to probe the LF 
over a wide redshift range and to sample the faint end of the 
LF down to unprecedented faint luminosities. We use photo¬ 
metric redshift spectroscopically calibrated with a sample of 
more than 230 spectroscopic redshift, ~ 80 of which in the 
HDF-S. We pay particular attention to the calculation of 
the rest-frame near-IR luminosities by comparing different 
methods and libraries of models. In §2 we describe the pho¬ 
tometric catalogue and the criteria adopted to construct the 
final sample of galaxies used to derive the LF. In §3 we de¬ 
scribe the procedure used to obtain calibrated photometric 
redshifts and we present the redshift and color distribution 
of galaxies. In §4 we discuss the derivation of the rest-frame 
Js-band and Ks-band absolute magnitudes and assess the 
influence of the different estimates on the LF. The resulting 
LF at various redshifts is derived in §5. In §6 our LFs are 
compared with the local estimates and with those previously 
obtained by other authors at comparable redshift to depict 
the evolution of the LF of galaxies up to 2 ~ 3. In §7 we 
summarize and discuss the results. 

Throughout this paper, magnitudes are in the Vega sys¬ 
tem unless explicitly stated otherwise. We adopt an = 
0.3, Ha = 0.7 cosmology with Ho = 70 Km s“^ Mpc“^. 


2 PHOTOMETRIC CATALOG AND SAMPLE 
SELECTION 

The Ks-band selected catalog we use has been obtained 
by combining the version 2 of the Hubble Space Telescope 
(HST) images in the U300, B 450 , Veoe and I 814 bands (Gaser- 
tano et al. 2000) with the deep (~ 30 hrs per filter) VLT- 
ISAAC observations centered on the HDF-S carried out in 
the framework of the FIRES project (Eranx et al. 2000; 
Labbe et al. 2003a) in the Js, H and Ks bands. The for¬ 
mal limiting magnitudes (3a within 0.2 arcsec^) reached are 
26.7, 29.1 29.5, 28.6, 26.2, 24.5 and 24.4 in the U300, B 450 , 
Veoe, I 8 i 4 , Js, H and Ks bands respectively. The detection 
has been performed using SExtractor (Bertin & Arnouts 
1996) on the Ks-band image. The Ks-band magnitude is the 
MAGJ3EST measured by SExtractor while the colors have 
been measured within the Ks-band detection isophote by 
using SExtractor in double image mode. The near-IR data 
reduction, the detection and the magnitude estimate are de¬ 
scribed in detail in Saracco et al. (2001; 2004). To construct 
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Figure 1. Upper panel — Ks-band number counts as resulting 
from the whole sample of 610 galaxies detected in the HDF-S. 
The counts raise till Ks~ 24 and drop at fainter magnitudes in 
agreement with the 100% completenes estimated at Ks~ 23.2. 
Lower panel — Photometric errors as a function of Ks-band mag¬ 
nitude for the whole sample of galaxies in the HDF-S. At Ks~ 23 
the photometric error is still lower than 0.2 mag for most of the 
galaxies while is much larger at fainter magnitudes. The dotted 
line marks the limiting magnitude of the K23 selected sample. 


our catalog, we considered the area of the HDF-S fully cov¬ 
ered by all the images in the various bands resulting in an 
area of ~ 5.5 arcmin^. On this area we detected 610 sources 
at Ks< 26. In Fig. 1 the observed number counts are shown. 

From this Ks-band selected catalog we extracted a com¬ 
plete sample of galaxies to derive the LF. On the basis of 
the simulations described in Saracco et al. (2001), we esti¬ 
mated that our catalogue is 100% complete at Ks~ 23.2. 
This is consistent with the number counts shown in Fig. 
1 which raise till Ks~ 24. Down to Ks~ 23 the near-IR 
photometric errors are lower than 0.2 magnitudes while 
they rapidly increase at fainter magnitudes as shown in 
the lower panel of Fig. 1 making very uncertain the pho¬ 
tometric redshift estimate. Thus, according to these con¬ 
siderations, we constructed the sample by selecting all the 
sources (332) brighter than Ks=23 (K23 sample hereafter). 
We then identihed and removed the stars from the K23 
sample relying both on the SExtractor star/galaxy classi¬ 
fier and on the colors of the sources. In Fig. 2 the stel¬ 
lar index (CLASS.STAR) computed by SExtractor for the 
332 sources of the K23 sample is shown as a function of 
the Ks-band apparent magnitude. Eor magnitudes brighter 
than Ks~ 21 point like sources (CLASS_STAR~ 1) and ex¬ 
tended sources (CLASS_STAR~ 0) are well segregated while 
they tend to mix at fainter magnitudes. We defined stellar 
candidates those sources with magnitudes Ks< 21 having a 
CLASS_STAR> 0.95 both in the Ks-band image and in the 
I 8 i 4 -band image. We obtained a sample of 23 stellar candi¬ 
dates. Among them, 7 bright (Ks< 18.5) sources display 
spikes in the HST-WFPC images implying that they are 


Figure 2. Sextractor stellar index (CLASS-STAR) as a function 
of Ks magnitude for the 332 sources brighter than Ks=23. 

most likely stars. For the remaining 16 candidates we ver¬ 
ified that their optical and near-IR colors were compatible 
with those of stars. In Fig. 3 the colors Js-Ks and Veoe- 
I8i4 of the K23 sample are plotted as a function of the Ks 
magnitude and of the Js-Ks color respectively. The 16 star 
candidates and the 7 bright stars are marked by filled points 
and starred symbols respectively. As can be seen, all but 2 
of the candidates occupy the stellar locus at J-K< 0.9 and 
are well segregated both in the color-magnitude and in the 
color-color plane. On the contrary, the two candidates ly¬ 
ing far from the stellar locus of Fig. 3 display colors not 
compatible with those of stars. In particular, while their J- 
K color (J-K> 1.3) could be consistent with the near-IR 
color of an M6-M8 spectral type, the other optical colors 
differ even more than 2 magnitudes from those characteriz¬ 
ing these stars. Thus, these two sources are most probably 
misclassified compact galaxies. This has made us to classify 
21 stars out of the 23 candidates originally selected on the 
basis of the SExtractor stellar index at Ks< 21. 

At magnitudes fainter than Ks=21 other sources (18) lie 
in the locus occupied by stars defined by a color J-K< 0.9. 
By a visual inspection of the 18 sources on the HST im¬ 
ages we clearly identified the extended profile for 7 of them. 
On the contrary, the remaining 11 sources have FWHM and 
surface brightness profile consistent with a point-like source. 
We probed the stellar nature of these 11 star candidates by 
comparing their observed spectral energy distribution (SED) 
defined by broad-band photometry with a set of spectral 
templates of stellar atmospheres from the Kurucz’s atlas 
(Kurucz 1993^). The observed SEDs of all the 11 star can¬ 
didates are very well fitted by those of M and K stars. All 
but one are not targets of spectroscopic observations, as we 
verified through the ESO archive data available. The spec¬ 
trum available for the remaining one (J32m58.07s-32’58.9’) 

1 www.stsci.edu/hst/observatory/cdbs/k93models.html 
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Figure 3. Colors Js-Ks (left) and V606-I814 (right) of the 332 Ks^ 23 as a function of the Ks magnitude and of the Js-Ks color 
respectively. The 25 star candidates are marked as filled symbols. Among them, we marked with starred symbols those sources clearly 
identified as stars because showing spikes in the HDF-WFPC images. 


confirms its stellar nature as M star. In the left panel of 
Fig. 4 we show the reduced of the best fitting template 
to the SED of the 18 sources with J-K< 0.9 obtained with 
the stellar atmosphere templates of Kurucz (Xstar) and with 
the stellar population templates of BC03 (Xgoi)- If is inter¬ 
esting to note the effectiveness of this method in identifying 
stars and galaxies among faint K-selected objects. In the 
right panel of Fig. 4 we show, as example, the SEDs of 5 
out of the 11 stars together with the relevant best-fitting 
Kurucz’s templates. Therefore, we classified stars also these 
11 sources and we removed a total of 32 stars from the K23 
sample which results in 300 galaxies. 


3 PHOTOMETRIC REDSHIFTS 
3.1 Templates selection 

Photometric redshifts have been derived by comparing the 
observed flux densities to a set of synthetic templates based 
on the latest version of the Bruzual & Chariot (2003) mod¬ 
els. The x^ minimization procedure of Hyperz (Bolzonella et 
al. 2000) has been used to obtain the best fitting templates. 
The set of templates used has been selected among a large 
grid of models to provide us with the most accurate estimate 
of the redshift for the galaxies in our photometric sample. 
The grid from which we selected the final set of templates, 
includes declining SFRs with time scale r in the range 0.1- 
15 Gyr besides the Simple Stellar Population (SSP) and the 
Constant Star Formation (cst) models. The templates have 
been produced at solar metallicity with Salpeter, Scalo and 
Miller-Scalo initial mass functions (IMFs). Various ranges 
of extinction Ay and two different extinction laws (Calzetti 
et al. 2000 and Prevot et al. 1984) have been considered in 
the best-fitting procedure. 

The selection of the best set of templates has been per¬ 


Table 1. Parameters defining the best set of templates obtained 
by comparing the photometric redshift with the spectroscopic red- 
shift of 232 galaxies, 151 in the HDF-N and 81 in the HDF-S. 


Best set 

SFHs T [Gyr] 

0.1, 0.3, 1, 3, 15, cst 

IMF 

Miller-Scalo 

Metallicity 

Z0 

Extinction law 

Calzetti 

Extinction 

0 < E{B - V) < 0.3 


formed by comparing the photometric redshifts Zp provided 
by the various sets of templates obtainable by the whole 
grid of models, with the spectroscopic redshifts Za of a con¬ 
trol sample of 232 galaxies. This sample includes 151 galax¬ 
ies in the HDF-N (Cohen et al. 2000; Dawson et al. 2001; 
Fernandez-Soto et al. 2002) and 81 in the HDF-S (Vanzella 
et al. 2002; Sawicki & Mallen-Ornelas 2003; Labbe et al. 
2003b; Trujillo et al. 2004; Rigopoulou et al. 2005) with 
known spectroscopic redshift. The best set of templates is 
the one minimizing the mean and the standard deviation 
of the residuals defined as Az = {zp — Zs)/p + Za) and the 
number of outliers defined as those sources with Az > 4 (ja^ . 
Both mean and standard deviation are iteratively computed 
applying a 4 ct clipping till the convergence. The final number 
of outliers includes all the galaxies removed by the clipping 
procedure. 

The set of templates and the range of extinction pro¬ 
viding the best results are summarized in Tab. 1. It is worth 
noting that the resulting best set does not include the SSP, 
contrary to some of the previous photometric redshift analy¬ 
sis. This is due to the fact that SSP templates always provide 
a better fit to the photometry then other SFHs but a wrong 
redshift estimate, typically lower than the real one. This re¬ 
sults both in a larger deviation of the residuals from the null 
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Figure 4. Left - Reduced obtained by fitting the SED of the 18 sources with K> 21 and J-K< 0.9 with the Kurucz’s templates (xltar) 
and with stellar population templates (x^^i)- The filled symbols mark the 7 galaxies. Right - The templates of stellar atmospheres from 
the Kurucz’s atlas are superimposed on the observed SEDs of 5 out of the 11 star candidates fainter than Ks> 21. 



Figure 5. - Left: Comparison of photometric and spectroscopic redshifts for 232 galaxies (upper panel), 151 of which in the HDF-N 
(empty circles) and 81 in the HDF-S (filled circles). The outliers are marked by dotted circles. The numbers refer to the ID number in 
the Fernandez-Soto et al. photometric catalog for the HDF-N and to our photometric catalog for the HDF-S. In the lower panel the 
residuals Az are plotted as a function of the spectroscopic redshifts. Right Distribution of the residuals Az (histogram) for the whole 
spectroscopic sample. A Gaussian (dashed line) with a = 0.065 and a mean deviation of (A^) = —0.015 from zero is superimposed on 
the observed distribution. 


value and in a larger scatter. In Fig. 5 the photometric red¬ 
shifts obtained with the best set of templates are compared 
to the spectroscopic redshift of the 232 galaxies of our con¬ 
trol sample. The typical scatter in our photometric redshift 
estimate is caz = 0.065 with a negligible deviation of the 
residuals from the null value ((Az) = —0.015). It can be seen 
from Fig. 5 that no correlation is present between the resid¬ 


uals and the redshift. We obtain three outliers (1%) over the 
redshift range 0 < z < 6, 2 out of which (IDs 892 and 242) 
in the HDF-S. Both these are well above the background on 
the I 8 i 4 -band image while their surface brightness appear 
very faint in the Ks-band image. This fact could have im¬ 
plied a wrong estimate of the colors and, consequently, of the 
photometric redshift. However, looking at the spectrum of 
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Figure 6. Cumulative distributions of photometric (continuous 
line) and spectroscopic (dashed line) redshift of the 232 galaxies 
with known redshift. The K-S test performed shows that the two 
distributions belong to the same parent population. 

892, we verified that its spectroscopic redshift is rather un¬ 
certain since it is based on a single emission line tentatively 
identified as Ha. It should be noted that both of them are 
fainter than Ks=23, the magnitnde at which we selected the 
complete sample to derive the LF (see §4 and §5). In Tab. 
2, we report the valnes of (Aa) and of caz together with 
the number of outliers relevant to the spectroscopic samples 
in the HDF-N, in the HDF-S and in the whole sample for 
different redshift bins. 

We verified the reliability of our photometric redshifts 
by comparing the spectroscopic and the photometric redshift 
distributions of the 232 galaxies. In Fig. 6 the cumulative 
distributions of the photometric and the spectroscopic red¬ 
shift are shown. A Kolmogorov-Smirnov test has provided a 
probability P(ks)— 40% that the two distributions belong 
to the same parent population confirming the reliability of 
our estimate. The same result has been obtained considering 
separately the 151 galaxies in the HDF-N and the 81 in the 
HDF-S. 

3.2 Redshift and color distributions 

Once defined the set of templates and parameters providing 
the most accurate estimate of the photometric redshift we 
have run the minimization procedure on the complete 
sample of 300 galaxies at Kss:23. For those galaxies having 
a spectroscopic redshift (74 at KsS(23), the best fitting has 
been constrained to the observed value. In Figure 7 the red¬ 
shift distribution of the sample is shown. The distribution 
has a median redshift Zm = 1-15 and has a tail extending 
up to Zj, ~ 6. There are 14 galaxies with photometric red¬ 
shift > 4, 5 of which at 5 < 2 :p ^ 6. The Js-Ks color of 
the 300 galaxies is shown in Fig. 8. Three of the Js-Ks>3 
galaxies were noticed on shallower near-IR images (Saracco 
et al. 2001) and previously analysed resulting in high-mass 



z 


Figure 7. Redshift distribution of the complete sample of 300 
galaxies brighter than Ks=23 in the HDF-S. The median redshift 
is 2 p = 1.15. The dotted line at Zp = 4 marks the redshift limit 
of the highest redshift bin considered in the LF estimate (see §5). 

evolved galaxies at 2 < z < 3 (Saracco et al. 2004). In Fig. 
8 we also plot the expected Js-Ks color obtained in the case 
of E(B-V)=0 assuming a declining star formation rate with 
time scale t = 1 Gyr seen at 5, 2 and 0.5 Gyr (from the top, 
thick lines) and a constant star formation (cst) seen at 1 
Gyr (dot-dashed line). The 1 Gyr old cst model is reddened 
by E(B-V)=0.5 to account for dusty star-forming galaxies 
(dot-dashed line). Most of the galaxies lie within the region 
limited by the 5 Gyr old model toward the red and by the 
0.5 Gyr old model toward the blue. The reddest galaxies of 
the sample (at 2 < z < 3) can be described by old stellar 
populations, i.e. populations with age comparable to the age 
of the Universe at the relevant redshift. 


4 K-CORRECTIONS AND NEAR-IR 
LUMINOSITIES 

The estimate of the k-correction and of the absolute magni¬ 
tude in the near-IR rest-frame is a critical step when deal¬ 
ing with high redshift galaxies. Indeed, unless having pho¬ 
tometry sampling the near-IR rest-frame, it implies an ex¬ 
trapolation with respect to the observed wavelength. This 
extrapolation is usually based on the best-fitting template 
which should be a good approximation of the true (but un¬ 
known) SED of the galaxy. In fact, the larger the redshift, 
the wider the extrapolation and, possibly, the uncertainty 
affecting the rest-frame near-IR luminosity. Moreover, this 
estimate can also be affected by systematics due to the dif¬ 
ferent libraries of models which can differ substantially in 
the near-IR domain. We have tried to assess whether and 
how such uncertainties affect the estimate of the rest-frame 
near-IR luminosity. We first checked the robustness of the 
rest-frame near-IR luminosity with respect to the extrap¬ 
olated photometry and than with respect to the different 

















The evolution of the near-IR LF of galaxies to z ^ 3 


7 


Table 2. Mean and standard deviation of the residuals Az for the spectroscopic samples of 232 galaxies (151 galaxies in the HDF-N 
and 81 in the HDF-S). For each estimate the results obtained without and with the 4cr clipping are shown. 



0 < 2 
Az ± (TAz 

< 2 

outliers 

2 < 2 
A 2 ± (TAz 

< 6 

outliers 

0 < 
A 2 ± (TAz 

z < 6 

outliers 

HDF-N 

0.0 ± 0.066 
0.0 ± 0.066 

0/118 


-0.044 ± 0.068 
-0.044 ± 0.068 

1/33 (3%) 

-0.009 ± 0.068 
-0.007 ±0.064 

1/151 

«1%) 

HDF-S 

0.033 ± 0.378 
0.033 ± 0.063 

2/71 

(3%) 

-0.017 ±0.075 
-0.017 ±0.075 

0/10 

0.026 ± 0.353 
0.030 ± 0.065 

2/81 

(2%) 

HDFs 

0.012 ± 0.230 
-0.011 ± 0.066 

2/189 

(1%) 

-0.037 ±0.070 
-0.037 ±0.070 

0/43 

0.002 ±0.210 
-0.015 ±0.065 

3/232 

(1%) 



Figure 8. Js-Ks color of galaxies versus redshift. The thick curves 
represent (from the top) the Js-Ks color derived by a declining 
star formation rate (e”*/'^) with e-folding time r = 1 Gyr and 
E(B-V)=0 seen at 5, 2 and 0.5 Gyr. The thin curve is a dusty 
(Av^=2) constant star formation (cst) seen at 1 Gyr. The tem¬ 
plates have been obtained with a Salpeter IMF at solar metal- 
licity. The dotted line at 2 : = 4 marks the redshift limit of the 
highest redshift bin considered in the LF estimate (see §5). 

models. For each galaxy, we compared the rest-frame Js- 
band absolute magnitude obtained from the equation 

MjRkj) = Js — 5log{DL{z)) - 25 - kj{z) (1) 

where Dl{z) is the luminosity distance at redshift 2 ; in units 
of Mpc, to the one derived by the equation 

MjRkjy) = my - 5log{DL{z)) - 25 - kjy{z) (2) 

In eq. 1 kj{z) = {Js,rest — Jsiz))temp is the conventional k- 
correction obtained as the difference between the rest-frame 
and observed magnitude computed on the best-fitting tem¬ 
plate. To this end, the transmission curve of the ISAAC- 
Js filter has been multiplied with the best fitting tem¬ 
plate at rest and redshifted to 2 (where z = Zs for the 74 
galaxies with known redshift in the K23 sample). In eq. 2, 
kjy = {Ja,rest — Y{z))temp and TUy is the observed apparent 


magnitude in the filter Y that best matches the rest-frame 
Js-band of the galaxy at the relevant redshift (e.g. Lilly et 
al. 1995; Pozzetti et al. 2003). In our case Y=Js, H, Ks 
according to the redshift of the galaxy. Both the “color k- 
correction” terms,fcj)^) and kjy, include the flux density 
dimming factor (1 -|- z) independent of wavelength. In Fig. 
9 (left panel) the difference between the rest-frame Js-band 
absolute magnitude derived for each galaxy from equations 
1 and 2 is plotted as a function of 2 . Up to 2 ~ 2, where 
the near-IR filters sample wavelengths Ajs > 0.4 gm and 
Aifs > 0.7 gm, the two equations provide the same values 
within an rms of about 0.15 mag, comparable to or lower 
than the typical photometric error. This dispersion does not 
affect the estimate of the LF. This is shown in the right panel 
of Fig. 9 where the distributions of the rest-frame Js-band 
absolute magnitude derived from eq. 1 (solid histogram) and 
from eq. 2 (dotted histogram) are compared for galaxies in 
different redshift bin. At 2 < 2, the two distributions be¬ 
long to the same parent population as confirmed by the K-S 
test we performed. At larger redshift, the conventional eq. 
1 provides absolute magnitudes which can differ systemat¬ 
ically from those derived from eq. 2. These differences can 
bias the estimate of the LF. The absolute magnitude dis¬ 
tributions shown in Fig. 9 for galaxies at 2 > 2 are indeed 
rather different, even if the difference is only marginally sig¬ 
nificant. Thus, until the near-IR band chosen to construct 
the LF (the Js filter in our example) samples the part of 
the SED red-wards A > 0.4 gm, characterized by a reg¬ 
ular shape, dominated by the emission of older stars and 
weakly affected by dust extinction and star formation, eq. 1 
and 2 provide consistent absolute magnitudes. In this case, 
having photometry approaching the rest-frame wavelength 
of the chosen filter (the Ks-band in our example) and us¬ 
ing eq. 2 do not improve the accuracy and the reliability of 
the estimate of the LF. On the contrary, when the chosen 
near-IR filter samples the blue and UV part of the spectrum 
(Ajs < 0.4 gm at 2 > 2), the observed emission in that filter 
(Js) is dominated by young (< 0.5 Gyr) stars which do not 
contribute significantly to the red part of the spectrum, i.e. 
to the rest-frame luminosity we want to derive. Moreover, it 
can be strongly affected by dust extinction. In this case the 
extrapolation to the near-IR rest-frame wavelength required 
by eq. 1 can give wrong values and the photometry in the 
redder band, the Ks filter in our case, is needed to apply eq. 
2 . 
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Figure 9. Left. The difference between the rest-frame Js-band absolute magnitude derived for each galaxy from eqs.l and 2 is plotted as 
a function of z. Right The distribution of the Js-band absolute magnitude of galaxies derived from eq. 1 (solid histogram) are compared 
to the one derived from eq. 2 (dotted histogram) for different redshift bin. No significant differences are present between the distributions 
at 2 < 2 as confirmed by the high probability provided by the K-S test (small sub-panels). 


We then checked the dependence of the LF on the dif¬ 
ferent library of models used to derive the best-fitting tem¬ 
plate. Recently, Maraston (2005; M05 hereafter) has shown 
the importance of the TP-AGB phase in modeling young 
stellar populations. This evolutionary phase can play an im¬ 
portant role in the continuum emission of a stellar popu¬ 
lation at A > 0.4 — 0.5 /rm with ages younger than 1 Gyr. 
The modeling of this phase in the Maraston’s code produces 
a brightening in the near-IR with respect to the same stel¬ 
lar populations (equal age, IMF and metallicity) modeled by 
the Bruzual & Chariot (2003) and PEGASE (Fioc & Rocca- 
Volmerange 1997) codes. To assess whether this difference 
affect the estimate of the LF, we compared the distributions 
of the Ks-band absolute magnitudes obtained with the BC03 
models and with the M05 models. We did not include the 
PEGASE models in this comparison since they provide spec¬ 
tral energy distributions very similar to those of BC03. We 
constructed the same set of templates described in Tab. 1 
using the SSPs of M05 and we searched for, for each galaxy, 
the best-fitting template at the relevant redshift. In Fig. 10 
the absolute magnitude distributions obtained with the two 
sets of models for the whole sample of galaxies at 2 < 4 and 
in three redshift bin are shown. The K-S test we performed 
does not show evidence of difference between the distribu¬ 
tions providing probabilities P(xs)> 0.8 in all the redshift 
bin. 

The above results show that our study of the LF and 
of its evolution with redshift is not dependent either on the 
method used to derive the rest-frame absolute magnitudes 
or on the library of templates used at least down to 2 ~ 2. 
Therefore, we are confident of the reliability of our estimate 
of the LF at least down to this redshift. However, we chose 
to push our study of the LF at 2 < 4 since the Ks-band still 


samples the red part {Xks > 0.45 fim) of the SED down to 
this redshift limit. We will use eq. 2 to derive the rest-frame 
Js-band absolute magnitudes and the conventional form of 
eq. 1 to compute the Ks-band absolute magnitudes. In Fig. 
11 the Ks-band k-correction derived for the sample of 285 
galaxies at 2 < 4 is shown as a function of redshift. For com¬ 
parison, superimposed on the data, we plot the k-correction 
as derived by 3 synthetic templates (continuous lines) ob¬ 
tained with a declining star formation rate seen at three 
different ages, (from the top) 5, 2 and 0.5 Gyr. Besides the 
models, we also plot (dashed line) the k-correction derived 
by the mean observed spectrum of local elliptical galaxies 
(Mannucci et al. 2001). 


5 THE NEAR-IR LUMINOSITY FUNCTION 
IN THE HDF-S 

5.1 Computing the LF with 1/Vmax method 

We computed the luminosity function, ${M)dM, using 
the non-parametric l/Vmax method (Schmidt 1968). This 
method does not require any assumption on the form of the 
LF and provides an unbiased estimate if the sample is not 
affected by strong density inhomogeneities. Previous studies 
of the LF of galaxies in the HDFs obtained with different 
estimators, both parametric and non-parametric (Takeuchi 
et al. 2000; Bolzonella et al. 2002), show no differences with 
respect to the LF obtained with the Vmax formalism. This 
implies that no strong inhomogeneities are present in the 
HDF-S and that our results are not dependent on the LF 
estimator used. In the 1/Vmax method, given a redshift bin 
[zi,Zu), each galaxy contributes to the number density of 
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Figure 10. The K-band absolute magnitudes obtained with the 
BC03 (Bruzual & Chariot 2003) models (solid histogram) for the 
sample of 285 galaxies at 2 : < 4 are compared with those obtained 
with the M05 (Maraston 2005) models (dotted histogram) for dif¬ 
ferent redshift bins. The distributions belong to the same parent 
population as shown by the high probability provided by the K-S 
test we performed (small panels). 
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Figure 11. Ks-band k-correction as a function of redshift for the 
285 galaxies brighter than Ks=23 at 2 ^ 4 in the HDF-S (small 
open circles). The filled squares are the median values of the k- 
correction in redshift bins 0.5 width (horizontal errorbars). The 
errorbars along the y-axis represent the scatter within each bin. 
The continuous lines are the k-corrections derived by a model 
described by an exponentially declining star formation rate with 
r = IGyr and t (from top) 5, 2 and 0.5 Gyr. The dashed line 
represent the k-correction of local ellipticals derived by the mean 
observed spectrum of ellipticals (Mannucci et al. 2001). 


Figure 12. Rest-frame Js-band luminosity function of galaxies 
derived with the 1/Vmax method (circles) in the three redshift bin 
considered: 0 < 2 < 0.8 {zm = 0.55, upper panel), 0.8 ^ 2 < 1.9 
(^Zm = 1-2, middle panel) and 1.9 ^ 2 < 4 {zm = 3.1, lower panel). 
We marked with open circles the LF obtained with the Js-band 
luminosities derived from eq. 1. The starred symbols show the 
LF uncorrected for the Ymax- The faint end (Mjg > —17) of the 
LF is defined by 9 galaxies at 2 < 0.3. Superimposed on the LF 
is shown the formal fit obtained with a Schechter function (thin 
line). The dotted line is the local LF of Cole et al. (2001). 


galaxies in that bin an amount inversely proportional to the 
maximum volume: 


K 


max.z • 


min{zu 


i) 


(fV 

dzdQ. 


dzdO. 


( 3 ) 


where Zmax,i is the maximum redshift at which the galaxy 
i with absolute magnitude Mi £ [M, M + dM] is still de¬ 
tectable given the limiting apparent magnitude Ks=23 of the 
sample, d^VjdzdO, is the comoving volume element and Q is 
the solid angle covered by the surveyed area. In our case, the 
sample has been selected over an area of about 5.5 arcmin^ 
of the HDF-S corresponding to a solid angle Q ~ 4.65 x 10“^ 
std. The comoving number density of galaxies ^{M)dM for 
each absolute magnitude bin in a given redshift bin and its 
error cr,j are computed as 


^{M)dM = V —2 —dM ai. 

' ^ ymax,i 

i 



We defined the different redshift bins so that a comparable 
and statistically significant number of galaxies fall in each of 
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Figure 13. Rest-frame Ks-band luminosity function of galaxies 
and Schechter best-fitting. Symbols are as in Fig. 12. 


them. According to this criterion and to the results obtained 
in §4, we divided the sample of 285 galaxies at 2 ; < 4 in 
the three redshift bins [0;0.8), [0.8;1.9) and [1.9;4). Each 
bin includes 101, 100 and 84 galaxies of which 52 (50%), 12 
(12%) and 10 (12%) with spectroscopic redshift respectively. 
The median redshifts in the three bin are Zm = 0.55, Zm = 
1.2 and Zm = 3.1 respectively. 

The rest-frame Js-band and Ks-band LFs we derived in 
each redshift bin are shown in Fig. 12 and 13 respectively 
(filled points). The starred symbols denote the LF uncor¬ 
rected for the incompleteness and show where the 1/Vmaa: 
correction takes place. In Fig. 12, the Js-band LF obtained 
with the rest-frame Js-band absolute magnitudes derived 
from eq. 1 is also shown (empty points). The figures suggest 
a systematic decrease of the number density of bright galax¬ 
ies coupled with a systematic brightening of the LF going 
to high redshift. This trend is present both in the Js-band 
LF and in the Ks-band LF. We will probe further these 
features in the next section. Besides this, the figures are 
also suggestive of a raise of the LF at faint absolute magni¬ 
tudes (Mjs > —17 and Mxs > —18) in the lowest redshift 
bin, raise which cannot be probed in the two bins at higher 
redshift where brighter galaxies are selectively missed. The 
same slope of the faint end is obtained by computing the LF 
in the redshift range 0 < 2 < 0.4 suggesting that the raise 
is not due to possible exceeding values of 1/Vmax- How¬ 
ever, it is worth noting that the raise we see is based on 
a very low statistics. Indeed, the faint end is defined by 9 
out of the 12 galaxies at 0.1 < 2 < 0.3 in our sample. It 


is indeed expected that, given the criteria adopted to select 
the HDF-S (devoid of bright sources), galaxies at low red¬ 
shift in this field are low luminosity galaxies (Mic >-20). 
This, coupled with the very faint limit (Ks=23) of the sam¬ 
ple, allows to probe the faint end of the luminosity func¬ 
tion at ( 2 ) ~ 0.2 down to the unprecedented faint near- 
IR magnitude M~ —14. These low luminosity galaxies are 
slightly bluer in the optical (0.45<V606-l8i4 < 0.9) than the 
bulk of the galaxies in the HDF-S while their optical-IR col¬ 
ors (1.4<l8i4-Ks< 2.7, 0.9<Js-Ks<1.5) are consistent with 
those of the bulk. All but one are brighter than V=25.8 and 
1=24.9 and are in the range 21<Ks< 22.9. Most of them 
show an irregular morphology while two of them appear 
very compact. In Fig. 14 the Veoe (upper panels) and the 
I 8 i 4 -band (lower panels) images of 7 out of the 9 galaxies 
at 2 < 0.3 fainter than M_ft:=-18 are shown. Both the ap¬ 
parent raise of the LF at faint luminosities and the irregular 
morphology of the galaxies populating the faint end are con¬ 
sistent with previous estimates of the local LF of galaxies 
at optical wavelengths. Such estimates show that the faint 
end of the near-IR LF is dominated by irregular and dwarf 
galaxies whose contribution to the comoving number den¬ 
sity increases going to lower luminosities as found for the 
local LF in the optical rest-frame (e.g. Marzke et al. 1994; 
Marzke et al. 1998; Zucca et al. 1997; Folkes et al. 1999). 
We estimated a comoving number density of galaxies with 
magnitude -19<Mic <-14 of Ufaint = 0.4 ± 0.27 Mpc“® at 
( 2 ) ~ 0.2. This value agrees with the one derived locally 
by Zucca et al. (1997) for galaxies M(,^. > —17 and it is 
nearly one order of magnitude higher than the comoving 
density of brighter galaxies {fibright = 0.039 ± 0.005 Mpc“® 
for Mx < —19). However, the very low statistic, the un¬ 
certainty in the photometric redshift estimate and the small 
volume sampled by the HDF-S, especially at low redshift, 
prevent us to strongly constrain the faint end of the LF. 
The possible raise we see at very faint near-IR luminosities 
needs a larger sample to be established. 

5.2 The evolution of the LF of galaxies in the 
HDF-S 

In Figs. 12 and 13 we plot the formal fitting to the LFs we 
obtained with a Schechter function (Schechter 1976). The 
fitting is performed over the whole range of absolute mag¬ 
nitude spanned by the LFs in each redshift bin. The results 
of the fitting are summarized in Tab. 3 where we report the 
values of the parameters a, M* and (j)* in each redshift bin. 
In Fig. 15 the errors contour at 68% confidence level for the 
joint parameters [a, M*] and [<()*, M*] of the Ks-band LF fit 
are shown. The decreasing number density and the possible 
brightening previously noticed in the LFs going to high red¬ 
shift are evident in the plot. However, given the statistical 
correlation between the Schechter parameters (a, M* and 
fp*), we will probe further on the significance of this evolu¬ 
tion by fixing the a parameter. Fig. 15 suggests also that 
we cannot constrain the faint end of the LF and thus the 
parameter a at high redshift. However, we tried to constrain 
the evolution of the LF shape by means of the K-S test (see 
e.g. Feulner et al. 2003) which uses all the absolute magni¬ 
tudes without binning. We derived in each redshift bin the 
cumulative LF (CLF) which gives the number of galaxies 
per unit volume brighter than M. In the upper panel of Fig. 
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Figure 14. Veoe-band (upper panels) and I 8 i 4 -band (lower panel) images of 7 out of the 9 galaxies at 2 < 0.3 fainter than =-19 
galaxies. The images are 3x3 arcsec centered on the source. The intensity of the images are optimized to show the galaxies. Galaxies are 
displayed from left to right on the basis of their Veoe-band apparent magnitude. 




Figure 15. Error contours at 68% confidence level for the parameters [o, M*] (left panel) and [</<*, M*] (right panel) of the Ks-band 
Schechter fit of Fig. 13. 


Table 3. Parameters a, M* and ip* of the Schechter function 
obtained by fitting the Js-band and the Ks-band LFs in the three 
redshift bin. 
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16 we show the three CLFs obtained in the Ks-band. Since 
the K-S test compares the CLFs normalized to unit, this test 
is not able to detect variations of (f>* ■ The normalized CLFs 


provide information on the relevant distribution of galaxies 
within the range of absolute magnitude considered, i.e. on 
the form of the LF. For this reason we expect that the K-S 
test is more sensitive to the variation of a which defines the 
form of the distribution rather than of M*. 

To quantify the sensitivity of the K-S test to the varia¬ 
tions of a and M* for samples of the same size of our sample 
we simulated and compared among them various samples of 
100 objects described by a Schechter LF with different values 
of a and M*. We have found that, at 90% confidence level, 
the minimum variations we are able to detect on a and M* 
on our samples are ±0.3 and ±0.8 mag respectively. In the 
lower panel of Fig. 16 the normalized CLFs obtained in the 
three redshift bin are compared two by two (lowest vs mid 
redshift bin (upper panel) and mid vs highest redshift bin 
(lower panel)). The K-S test does not point out significant 
differences among the distributions as shown by the proba¬ 
bilities Pxs —0.61 (bin 1 vs bin 2) and Pxs —0.24 (bin 2 
vs bin 3). Thus, we conclude that if a changes its variation 
has to be lower than 0.3 from Zm. — 0.6 to Zm — 3. 

We then probed the evolution of M* and tp* by compar- 
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Figure 16. Left - Cumulative luminosity function of galaxies derived with the llVmax in the three redshift bin: 0 < z < 0.8 (circles), 
0.8 ^ z < 1.9 (squares) and 1.9 ^ z < 4 (triangles). Right - Two by two comparison of the cumulative LFs in the three redshift bin 
normalized over the range of absolute magnitude in common. 




Figure 17. Left - Superimposed on the Js-band LFs are shown the relevant best fitting obtained with a Schechter function with a = —1.0. 
Right - Same as in the left panel but for the Ks-band. 


ing the values obtained by fitting a Schechter function for a 
fixed value of a. We assumed the value a = —1.0 since it is 
very close to the values we derived in each redshift bin and 
to the values previously derived by other authors for local 
galaxies (e.g. Kochanek et al. 2001; Loveday 2000). The re¬ 
sults of the fitting are summarized in Tab. 4 and shown in 
Fig. 17. The possible brightening combined to the decreasing 
number density of bright galaxies previously noticed in the 
LF from Zm — 0.6 to Zm — 1-2 is confirmed and, possibly, it 
extends to ~ 3. The normalization of the LF and, con¬ 
sequently, the number density of bright galaxies decreases 
by a factor 3 from Zm — 0.6 to Zm — 1.2. In parallel, the 
characteristic magnitude M* brightens by ^ 0.4 mag. The 
evolution of fp* is detected at high confidence level (> Scr), 
as shown by the error contours shown in Fig. 18 while the 


significance of the brightening of M* is Icr. The same trend 
persists, even if at a low confidence level for cp* (< 2a), from 
Zm — 1.2 to Zm — 3 where <p* {z ^ 3) ~ 0.5</)*(z ^ 1.2) 
and M* brightens by ~ 0.4 mag. It is worth noting that 
this evolution is not dependent on the value assumed for a. 
Indeed, by assuming a = —1.1 we obtain for the Ks-band 
LF fit M*(z ~ 0.6) = -24.62 and (p*{z - 0.6) = 0.0045 
Mpc“® to be compared with M*(z ^ 1.2) = —25.03 and 
<p*{z ^ 1.2) = 0.0016 Mpc“®, resulting in the same luminos¬ 
ity and density evolution (see Tab. 4). 












The evolution of the near-IR LF of galaxies to z ^ 3 13 


Table 4. Parameters M* and <j>* of the Schechter function fitting 
the Js-band and the Ks-band LF with a = —1.0 and the Ks-band 
LF with a = —1.1. to the LF in the three redshift bin. 
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Figure 18. Error contours at Icr and 3a for the parameters 0* 
and M* of the Schechter fit with a = —1.0 to the Ks-band LF. 

6 THE EVOLUTION OF THE NEAR-IR LF TO 

Z~ 3 

In this section we compare our results with those obtained 
by other authors both at lower and at comparable redshift in 
order to constrain the evolution of the near-IR LF of galaxies 
from 2 ^ 3 to 2 ~ 0. 


6.1 The Ks-band LF 

In Tab. 5 the parameters of the Schechter function of the Ks- 
band LF obtained by various authors with different samples 
are summarized. In Fig. 19 the fitting Schechter functions 
obtained by the various authors are shown in the range of 
absolute magnitudes reached by the surveys according to the 


values reported in Tab. 5 (Mx(min)). We first compared the 
Ks-band LF we derived at Zm — 0.6 with those derived on 
local samples of galaxies in order to probe and to constrain 
the evolution of the LF on the last 4 Gyrs. In particular, 
we considered the LFs of Cole et al. (2001), Kochanek et al. 
(2001) and that of Loveday (2000). These three local LFs 
are consistent among them. The former two are based on 
the largest samples of local galaxies while the latter extends 
down to very faint luminosities (Mx(min)~ —16). In Fig. 
19 the three local LFs are superimposed to the LF we de¬ 
rived at Zm — 0.6 (left upper panel) based on 101 galaxies 
50% of which with spectroscopic redshift. The agreement 
with the LFs of Kochanek et al. and of Loveday et al. is 
rather good. The largest deviation is respect to the LF of 
Cole et al. However, this deviation is not significant as also 
confirmed by the K-S test which gives Pks — 0.5. The char¬ 
acteristic magnitude of the Schechter function we derive at 
Zm ~ 0.6 are consistent with the local value (see Tab. 5) 
even is systematically brighter of about 0.2-0.3 magnitude 
while no differences are found in the normalization 0*. Thus, 
the comparison with the local near-IR LF of galaxies does 
not point out evidence of strong luminosity and/or density 
evolution at z < 0.8. Given the statistical errors of our esti¬ 
mate we conclude that the near-IR LF has evolved not more 
than 0.2-0.3 magnitude in the last 4 Gyr. 

We then compared our LF at Zm ~ 0.6 with those de¬ 
rived at comparable redshift. In the right upper panel of Fig. 
19 the Schechter functions representing the LFs obtained by 
Glazebrook et al. (1995), Feulner et al. (2003), Bolzonella 
et al. (2002) and Pozzetti et al. (2003) at 2 ^ 0.5 are su¬ 
perimposed to our estimate. Glazebrook et al. and Feulner 
et al. derive the LF assuming a value of a of -1.0 and -1.1 
respectively. The LF we obtained with these values of a, 
summarized in Tab. 4, are in very good agreement with the 
results of Glazebrook et al. and of Feulner et al. The largest 
deviations are respect to the LF of Bolzonella et al. and of 
Pozzetti et al. Both these LFs are characterized by values 
of a steeper and by brighter characteristic magnitudes with 
respect both to our and to the other LFs at these and at 
lower redshift. The K-S test we performed shows that the 
deviation from the LF of Pozzetti et al. is not significant 
{Pks — 0.15) while it is significant (Pks — 0.02) the devia¬ 
tion from the LF of Bolzonella et al. This latter result, rather 
surprising since we deal with the same field, is also evident 
by comparing the LF we obtained at higher redshift. In the 
lower left panel of Fig. 19 the LFs obtained by Pozzetti et al. 
at 0.75 < 2 < 1.3, by Kashikawa et al. (2003) at 1 < 2 < 1.5 
and by Bolzonella et al. at 1 < 2 < 2 are superimposed to 
the one we obtained in the redshift range 0.8 Si 2 < 1.9. 
The agreement with the LF of Pozzetti et al. is remarkable. 
On the contrary, it is evident that our LF is not consistent 
with the LFs of Bolzonella et al. and of Kashikawa et al. 
who find a very steep {a < —1.35) LF coupled with a very 
bright (M)(- < —25.5) characteristic magnitude. In both the 
cases the K-S test gives a probability less than 1% that the 
distributions are drawn from the same parent population. 
This is shown in Fig. 20 where the cumulative distributions 
relevant to these LFs are shown. The reasons of this dis¬ 
agreement have to be searched both in the different samples 
and magnitudes used by Bolzonella et al. and Kashikawa et 
al. 

As to the LF of Bolzonella et al., even if relevant to the 
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Figure 19. The rest-frame Ks-band luminosity function derived with the llVmax (points) and LF Schechter fitting (solid curve) are 
compared with the LFs obtained by various authors (Koch.= Kochaneck; Love.= Loveday; Glaz.= Glazebrook; Feul.= Feulner; Bolz.= 
Bolzonella; Pozz.= Pozzetti; Kash.= Kashikawa). In the upper-left panel the LF we derived in the redshift bin Zm ~ 0.6 is compared 
with the local K-band LFs. In the other panels our estimates are compared with those obtained by other authors at comparable redshift: 
z ~ 0.5 (upper right), z cz 1 (lower left) and z ~ 3 (lower right). 


HDF-S, it is based on a near-IR sample extracted from an 
optically selected catalogue (Vanzella et al. 2001). At vari¬ 
ance with respect to our sample, detection and magnitude 
estimates were performed and optimized on optical HST im¬ 
ages. Moreover, the selection and the completeness of the 
near-IR sample extracted by Bolzonella et al. are based on 
the I 8 i 4 -Ks color distribution. The near-IR data are 4 times 
shallower than those used in the present work and, conse¬ 
quently, the sample is ~ 1 magnitude shallower. It should 
also be noted that 25% of our sample (50% of the sample at 
2 < 0.8) has spectroscopic redshift and that these redshifts 
have been used both to optimize the photometric redshift 
and to estimate the LF. All that implies differences in the 
near-IR magnitude estimates which can account for the dif¬ 
ferent LF obtained. 

As to the comparison with the LF of Kashikawa et al., 
no obvious reasons can be put for the disagreement we ob¬ 
tained in this redshift bin. Cosmic variance affecting such 
small areas can affect the two samples accounting for this 
discrepancy. In the highest redshift bin considered, Zm — 3, 


the disagreement is even larger as shown in the lower-right 
panel of Fig. 19. At this redshift, the Ks-band samples rest- 
frame wavelengths A ~ 0.55 /im and the large extrapolation 
needed to derive the rest-frame Ks-band luminosity can be 
strongly dependent on the best-fitting template. In this case, 
as suggested by the analysis we presented in §5, different 
templates could imply different LFs. However, it should be 
noted that the disagreement is based on the comparison of 
our LF with the Schechter fit found by Kashikawa et al. 
and not with their data-points. The few data-points defin¬ 
ing their LF in the two highest redshift bin (Fig. 5 of their 
paper) suggest that the Schechter parameters could not be 
so strongly constrained and that the disagreement could be 
less severe. 

The agreement with the Ks-band LFs derived by other 
authors at redshift comparable to z ~ 0.6 and to 2 ~ 1.2 
confirms the evolution we detect in the LF of galaxies 
in the HDF-S, i.e. a brightening AM* ~ —0.4 coupled 
with a decrease of the number density of bright galaxies 
A(j)*/(j>* ~ —0.65 from Zm ~ 0.6 to Zm ~ 1.2 for a = —1.0. 
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Figure 20. Comparison among the cumulative LFs derived at 
2 ~ 1.2 normalized over the same range of absolute magnitude. 
The solid histogram is the CLF we derived with the 1 jNmax. ■ The 
probability that the cumulative distributions derived by the LF 
of Bolzonella et al. (dotted line) and of Kashikawa et al. (dashed 
line) are drawn from the population our LF belongs to is less than 
1 %. 

At larger redshift, the uncertainties in the estimate of rest- 
frame Ks-band luminosities due to the extrapolation based 
on the best-fitting template could affect the LF estimates 
and mid-IR observations would be useful. However, even 
if our LF at 2 m — 3 deviates from the other LFs at this 
redshift, these latter LFs trace the same trend we observe 
suggesting that the evolution we detect at 2 < 2 extends to 
2 > 3. 

6.2 The Js-band LF 

Analogous results are obtained by the comparison of our 
Js-band LF estimates with those previously obtained at 
comparable redshift. The comparison is shown in Fig. 21 
and extends down to 2m ~ 1.2 since no Js-band LFs have 
been estimated at higher redshift by other authors. The lo¬ 
cal Js-band LF of Cole et al. (2001; M} = —23.13 and 
F = 0.4 • 10“^ Mpc“®) and the local LF of Balogh et 
al. (2001, M*j = —23.02) are consistent with our LF at 
2m ~ 0.6. Thus, as in the case of the Ks-band, we do not 
find evidence of strong evolution of the LF at 2 < 0.6. How¬ 
ever, in agreement with the Ks-band LF, we point out a 
brightening of of about 0.2-0.3 magnitude from 2 = 0 
to 2m ~ 0.6. 

The Js-band LFs obtained at 2 ~ 0.5 by the various 
authors show a much larger scatter than those obtained in 
the Ks-band at comparable redshift. As for the Ks-band, 
our estimate agrees with the LFs of Pozzetti et al. and of 
Feulner et al. (both based on spectroscopic redshift) while 
deviates significantly from the LFs of Bolzonella et al. and 
of Dahlen et al. (2005). In particular, both these LFs are 
characterized by significantly fainter M* (0.4 mag with re¬ 
spect to our LF and 0.5-0.6 mag with respect to the LF of 
Feulner et al. and Pozzetti et ah). Also at redshift 2m ~ 1.2 


our LF agrees very well with the Pozzetti et al. LF while 
deviates significantly from the LFs of Bolzonella et al. and 
Dahlen et al. These latter find also an evolution of M* in 
the opposite sense with respect to what we and other au¬ 
thors found, i.e. a dimming from 2 ~ 0.4 to 2 ~ 0.9. They 
use eq. 2 to derive the rest-frame absolute magnitudes. They 
suggest that a negative evolution should be expected in Js 
and at longer wavelengths by extrapolating the results ob¬ 
tained by Ilbert et al. (2005), who find that the evolution 
of M* is systematically weaker from U to I-bands. Thus, a 
turnover in the evolution of M* could be expected at ~J and 
a negative evolution should take place at longer wavelength. 
They suggest that the positive evolution which is instead 
observed by many authors is explained by the method used 
to derive the rest-frame luminosities we discussed in §4: the 
observed band (Ks or Js) samples shorter rest-frame wave¬ 
lengths (where the evolution is stronger) at higher redshift 
and the use of eq. 1 to derive the rest-frame luminosities 
could mimic a brightening in M*. However, we have shown 
in §4 and §5 that this is not true at least down to 2 ^ 2 and 
we can consequently exclude with certainty that the evolu¬ 
tion found by us and by other authors is due to a wrong 
derivation of the rest-frame luminosities. 

Thus, from the comparison with the other LFs in Js- 
band we confirm an evolution in the LF of galaxies from 
2m ~ 0.6 to 2m ~ 1.2, i.e. a brightening AM* ~ —0.4 
coupled with a decrease of the number density of bright 
galaxies. This evolutionary trend extends to 2 ~ 3 even if, 
as for the Ks-band LF, mid-IR observations would be needed 
to confirm the trend. 


7 SUMMARY AND CONCLUSIONS 

We have probed the evolution of the rest-frame Js-band and 
Ks-band LFs of field galaxies with a complete sample of 
about 300 galaxies selected in the HDF-S at KsS( 23 (Vega). 
Photometric redshifts have been obtained from template 
SED fitting to U 300 , B 450 , Veoe, I 8 i 4 , Js, H and Ks photom¬ 
etry after having calibrated and optimized the procedure 
and the set of templates with a control sample of 232 spec¬ 
troscopic redshifts, 151 from the HDF-N and 81 from the 
HDF-S. The accuracy in the redshift estimate we obtained 
is 0.06 (rms). 

We investigate the reliability of the rest-frame near- 
IR absolute magnitudes obtained using the conventional 
method based on the extrapolation of the observed photom¬ 
etry on the best-fitting template and using the photometry 
approaching the rest-frame near-IR. We find that the rest- 
frame Js-band absolute magnitudes obtained through the 
photometry in the redder bands (H and Ks according to the 
redshift of the galaxy) are consistent with those obtained 
from the Js-band photometry at least down to 2 ~ 2. This 
shows that the LF is not dependent either on the extrapo¬ 
lation made on the best-fitting template or on the library of 
models used. 

We derived the Js-band and the Ks-band LF in the 
three redshift bins [0,0.8), [0.8,1.9) and [1.9,4) centered at 
the median redshift Zm ~ 0.6, Zm ~ 1.2 and 2 m ~ 3. Each 
bin contains 101, 100 and 84 galaxies of which 52 (50%), 12 
(12%) and 10 (12%) with spectroscopic redshift respectively. 
The analysis of the observed LF at different redshifts and 
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Table 5. Summary of the parameters of the Schechter function obtained by various authors by fitting the K-band LF of different samples 
of field galaxies. For each sample the relevant limiting magnitude, the area, the range of redshift covered by the survey (or the mean 
redshift of the sample) and the number of objects are reported. The parameters of the Schechter function and the faintest absolute 
magnitude sampled by the data are scaled to the cosmology adopted here {Q,m = 0.3, Ha = 0.7, Hq = 70 Km s“^ Mpc“^). 


Source 

‘^lim 

Area 

(arcmin^) 

z 

obj.s 

M* 

o. 

P 

10“® [Mpc“®] 

Mat (min) 

Mobasher et al. 1993^ 

K<13 Bj ^17 


0 -0.1 

181 

-24.14 ±0.30 

-1.00 ±0.30 

3.8 ±0.6 

-22 

Glazebrook et al. 1995 

K^17.3 

552 

0 -0.2 

55 

-23.91 ±0.23 

-1.04 ±0.30 

7.6 ± 1.8 

-22 




0 -0.8 

124 

-24.42 ±0.11 

-1.00 

4.1 ±0.4 

-23.5 

Cowie et al. 1996 

K<20 

26.2 

0-1 

393 

-24.27 

-1.25 

1.0 

-21 

Gardner et al. 1997 

K^15 

15840 

(0.14) 

465 

-24.07 ±0.17 

-0.91 ±0.24 

4.9 ±0.7 

-21.5 

Szokoly et al. 1998 

K<16.5 

2160 

0-0.4 

175 

-24.57 ±0.30 

-1.30 ±0.20 

3.9 ± 1.0 

-21.5 

Loveday 2000^ 

bj <17.15 


(0.05) 

345 

-24.35 ±0.42 

-1.16 ±0.19 

4.1 ±0.3 

-16 

Kochanek et al. 2001 

K<11.25 

< 2.5 X 10'^ 

(0.02) 

3878 

-24.16 ±0.05 

-1.09 ±0.06 

4.0 ±0.3 

-21 

Cole et al. 2001 

K<13 

2.2 X 10® 

(0.05) 

17173 

-24.21 ±0.03 

-0.96 ±0.05 

3.7 ±0.6 

-20 

Feulner et al. 2003 

K<17.5 

649 

0.1-0.3 

157 

-24.56 ±0.24 

-1.10 

3.8 ±0.4 

-21 




0.3-0.6 

145 

-24.81 ±0.26 

-1.10 

2.4 ±0.9 

-23.5 

Huang et al. 2003 

K<;i5 

1056 

(0.14) 


-24.47 ±0.08 

-1.39 ±0.09 

4.5 ±0.7 

-20.5 

Pozzetti et al. 2003 

K<20 

52 

0.2-0.65 

132 

-24.87 ±0.63 

-1.25 ±0.22 

1.8 ± 1.2 

-19.5 




0.75-1.30 

170 

-24.77 ±0.49 

-0.98 ± 0.45 

2.9 ± 1.4 

-23 

Drory et al. 2003® 

K<19.5 

998 

0.4-1.2 

~ 5000 




-23.5 

Bolzonella et al. 2002^ 

K<22 


0-1 


-25.03 ±0.28 

-1.17 ±0.09 

3.4 ± 1.2 

-17 




1-2 


-25.71 ± 1.14 

-1.42 ±0.10 

0.8 ±0.7 

-20.5 

Caputi et al. 2005^ 

K<22 

50 

0-2.5 

1600 





Kashikawa et al. 2003^ 

K<24 


0 .6-1 

439 

-25.41 ±0.20 

-1.35 ±0.04 

2.3 ±0.7 

-18.5 




1-1.5 


-25.53 ±0.27 

-1.35 ±0.06 

1.9 ±0.8 

-19.5 




1.5-2.5 


-25.58 ±0.12 

-1.37 ±0.05 

1.3 ±0.5 

-19.5 




2.5-3.5 


-24.66 ±0.23 

-1.70 ±0.11 

1.0 ±0.7 

-20 

This work® 

K<23 

5.5 

0 -0.8 

101 

-24.47 ±0.63 

-1.14 ±0.11 

4.0 ±2.0 

-14 




0 .8-1.9 

100 

-24.53 ±0.37 

-0.85 ±0.17 

2.9 ± 1.0 

-19 




1.9-4 

84 

-25.18 ±0.53 

-1.10 ±0.31 

1.0 ±0.6 

-22 


^ Based on K-band imaging of an optically selected sample 

^ Based on K-band data of an optically selected sample and on photometric redshift. 

® Based on photometric redshift. The number of galaxies (439) is the total number of galaxies in the sample. 


the comparison with those previously found by other authors 
at comparable redshift and at 2 = 0 provided the following 
results: 

1. We find hints of a raise of the faint end (Mjs > —17 and 
Mxs > —18) of the near-IR LF at Zm ^ 0.6. The raise is 
defined by 9 galaxies at 2 < 0.3 with irregular morphology 
as found for the local LF at optical wavelength (e.g. Marzke 
et al. 1994; Marzke et al. 1998; Zucca et al. 1997; Folkes 
et al. 1999). They account for a comoving number density 
h = 0.4 ± 0.23 Mpc“® at < 2 >~ 0.2, almost one order of 
magnitude higher than that of brighter galaxies. However, 
given the low statistics, such raise of the faint end needs a 
larger sample to be established; 

2. We find no evidence of a steepening with redshift of the 
near-IR LFs of galaxies in the HDF-S. The value of a we find 
is consistent with the local value. Given the size of our sam¬ 
ple, we estimate that if a changes with redshift its evolution 
is constrained within ±0.3 from 2 ~ 0 to 2 ~ 3. 

3. We do not find evidence of strong evolution of the LF up to 
Zm — 0.6, where 50% of our galaxies has spectroscopic red¬ 
shift. The comparison with the local estimates of the near-IR 
LF (Cole et al. 2001; Kochanek et al. 2001) shows that M* 
has evolved not more than 0.2-0.3 magnitude in this red¬ 
shift range, while the number density of galaxies {4>*) has 
not evolved significantly in this redshift range, in agreement 
with the previous studies (e.g. Drory et al. 2003; Pozzetti et 


al. 2003; Feulner et al. 2003; Glazebrook et al. 1995; Caputi 
et al. 2004, 2005). 

4. We clearly detect an evolution of the LF at Zm > 0.6 char¬ 
acterized by a brightening of M* of about 0.6 and a decrease 
of <j)* by a factor 2-3 to Zm ^ 1.2 both in the Js-band and in 
the Ks-band. By fixing a — 1.0, this evolution is character¬ 
ized by AM* ~ —0.6 (at Icr) coupled with A(j>*/(p* ~ —0.65 
(at > 3 (t). The brightening persists (at Icr) up to 2 „, ~ 3 
together with the decline of <j)* (at ~ 2cr) even if at a lower 
extent. 

Our results agrees, at least qualitatively, with most of 
the analysis previously done (Drory et al. 2003; Pozzetti et 
al. 2003; Feulner et al. 2003; Caputi et al. 2004, 2005) while 
deviates from the recent estimate of Dahlen et al. (2005). 
Taking into account that the near-IR light is rather tightly 
connected to the stellar mass, our results can give insights on 
the evolution of the stellar mass of the galaxies besides their 
luminosity. Up to 2 ~ 0.8, clearly little luminosity evolution 
and no density evolution is observed. This suggests that the 
population of local bright galaxies was already formed at 
2 ^ 0.8 and that their stellar mass growth was already com¬ 
pleted at this redshift. This agrees with the results on the 
evolution of the stellar mass density obtained by many au¬ 
thors: Fontana et al. (2003) and Rudnick et al. (2003) in the 
HDF-S, Fontana et al. (2004) on the K20 sample (Cimatti 
et al. 2002), Yamada et al. (2005) in the Subaru Deep Sur- 
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Figure 21. The rest-frame Js-band LF derived with the 1/Ymax 
(points) and LF Schechter fitting (solid curve) are compared with 
the LFs obtained by various authors. In the upper panel the LF 
in the redshift bin Zm — 0.6 is compared with the local J-band 
LFs. In the other panels our estimates are compared with those 
obtained by other authors at comparable redshift: z ~ 0.5 (middle 
panel), z ~ 1 (lower panel) 


vey Field, Bundy et al. (2005) and Drory et al. (2005) in 
the GOODS and FDF fields and Feulner et al. (2005a) in 
the MUNICS fields. Indeed, all of them find very little or 
even no evolution in the stellar mass function of galaxies 
and in the specific star formation rate (SSFR) at 2 < 1 for 
high-mass galaxies. 

At 2 > 0.8 the evolution is stronger and the larger 
brightening observed is accompanied by a decrease of the 
number density of bright/massive galaxies. In the redshift 
bin [0.8,1.9) the number density of bright galaxies is ~ 
30 — 50% of the local value and reaches ~ 20 — 30% at 
1.9 < 2 < 4. This decline in the number density suggests 
that up to 70% of the local massive/bright galaxies has 
grown atl<2<2 — 3 through star formation, merging 
or both. Thus, we should expect to observe massive star¬ 
forming galaxies with disturbed/irregular morphology and 
merging systems in this high redshift range, as indeed found 
in same cases (e.g. Daddi et al. 2004). On the other hand, the 
observed evolution implies also that at least 30% of the local 
bright/massive galaxies was already in place at 1 < 2 < 2 as 
indeed found by some authors (Fontana et al. 2004; Glaze- 
brook et al. 2004; Forster Schreiber et al. 2004). Pozzetti et 
al. (2003) and Gaputi et al. (2005a) show that the bright end 
of the LF is dominated by early-type galaxies and that their 
contribution to the bright end does not decline with redshift. 


Thus, that 30% should be dominated by early-types imply¬ 
ing that we should observe massive evolved galaxies fully 
assembled at 1 < 2 < 2 and even at 2 > 2. This is indeed 
observed by many authors (Gimatti et al. 2004; McGarthy 
et al. 2004; Saracco et al. 2004, 2005; Longhetti et al. 2005; 
Daddi et al. 2005; Labbe et al. 2005). It is worth noting that 
this picture could be accounted for, at least qualitatively, in 
the hierarchical picture of galaxy formation as suggested by 
the recent results obtained by Nagamine et al. (2005). 

From these considerations we can gather that the 
growth of massive galaxies does not follow a unique way but 
displays different behaviours. A significant fraction (50-70%) 
of the brighetst/most massive galaxies increases their stel¬ 
lar mass over a large redshift range at 2 > 1. The remaining 
fraction reaches their final mass in a narrower redshift range 
at 2 > 3 since it is already in place by this redshift. This 
suggests that, for the former, the stellar mass growth has 
been less efficient and has proceeded more slowly than for 
the latter. On the contrary, in the latter, the stellar mass has 
to be grown rapidly in a short interval, surely much shorter 
than 1 Gyr given their redshift. This is supported by the 
recent results derived by mid-IR massive galaxies observed 
at high-z (e.g. Gaputi et al. 2005b) and suggests a high effi¬ 
ciency in the accretion of the stellar mass in massive haloes 
in the early Universe, possibly through a very efficient star 
formation. The recent results obtained on the evolution of 
the SSFR of galaxies seems to point toward this direction. 
The strong increase in the SSFR of the most massive galaxies 
with redshift (Feulner et al. 2005b) favors indeed an efficient 
SFR in the brightest galaxies at 2 > 3 — 4 constraining their 
growth in a very short interval. 
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